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EXACT SOLUTION OF SOME LINEAR MATRIX EQUATIONS USING ALGEBRAIC METHODS 

Abstract 


Let PA+BP * -C be a linear matrix equation where A is an n*n matrix, 
B an m*m matrix and C an m*n matrix, all three matrices taken over the 
reals. Let R[x,y] be the ring of polynomials in two indeterminants x and 
over the reals R, and MW the set of all m*n matrices over the reals. Let 
4* * (<t> 2 (x), ^ 2 be tiie ideal in R tx,yl generated by 2 (x) and (y) , 
where <t> 2 (x) is the characteristic polynomial of A and the charac- 

teristic polynomial of B. The elements of the quotient ring R[x,y]/'F are 
denoted by 4' + a(x,y) * {h(x,y) + a(x,y) |h(x,y) e 4'}. 

Define the action f_.:Rtx,y] x MW *♦* MW in the following manner: 

BA 

f BA <h(x,y),M) = I h B j .M*A k 
jk 3 

where h(x,y) * £ h x^y^ 811 element in R[x,yJ , 
jk 3K 

M is an element in MW 


The action f allows for the interpretation of MW as a module over 
Rlx/yl/y with multiplication {*) of elements in R^yl/V with elements in 
MW given by: 

(4? + a(x,y) ) *M * f (a(x,y)mod 4», M) 

dA 

where is the element of minimal degree in ¥ + a(x,y) . 

The polynomial x+y is in the ooset f + (x+y) . In the event that 
4' + (x+y) has an "inverse" 4* + q u (x,y) in Rtxjyl/’F such that (f + (x+y))* 
(4* + q u (x,y)) ■ 4' + u where ¥ + u is the ooset containing the real number 


u f 0 we can write: 
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(¥ + (x+y))*P ** PA + BP = -C 

u-P - (f + q u (x # y))*(-C) * ^(^(XjyJmod ¥,-0 
P » “ • f BA (q u (X ' y,m ° d _C ** 
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J * 1 . Introduction 

In the past fifteen years or so there has been impressive progress in 
the theoretical understanding of the structure, representation and control 
of linear multivariable systems. In contrast, workers in the field have 
paid little attention to the computational aspects of systems problems. 
This does not mean that algorithms for the solution of systems problems 
have not been developed. But most of the algorithms that have been pro- 
posed have never been seriously studied as far as stability convergence 
and similar issues are concerned. Even the LQG problem, bulwark of the 
so called "modern control theory" seems to be little understood from the 
computational point of view. 

In this paper we undertake a study of solution methods for Linear 
Matrix Equations including Lyapunov's equation 

PA + A'P - -Q (1.1) 

using methods of modern algebra. The emphasis is on the use of finite 
algebraic procedures which are easily implemented on a digital computer 
and which lead to an explicit solution to the problem. 

It is well known that the Lyapunov equation is important in the study 
of stability of linear finite dimensional time-invariant systems. If £ .s 
symmetric and positive definite and if A is a stability matrix then the 
unique solution to (1.1) is given by the convergent integral 

r A't At 

P - I e Qe dt 

0 

(Cf . B ROCKETT) . 

However, the solution requires the evaluation of an integral over an 


' t t 1 t I I 1 
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infinite time interval. Other methods of solution do exist all of which 
have the basic drawback of giving an approximate solution. This becomes 
frustrating when the problem is ill-conditioned. 

The need for solving the Lyapunov equation also arises when one uses 
Newton's Method to solve the Algebraic Riccati equation (cf. KLEINMAN) 

(A - BR* 1 B'P)'P + P (A - BR -1 B'P) - -C'C - PBR _1 B'P. 

Here a finite number of Lyapunov equations have to be solved. 

This paper has been inspired by an important paper by KALMAN. Kalman's 
concern was the characterization of polynomials whose zeroes lie in certain 
algebraic domains (and the unification of the ideas of Hermite are Lyapunov) . 
In this paper we show that the same ideas lead to finite algorithms for 
the solution of linear matrix equations. 

This paper is divided into five sections. In section 2 we introduce 
the action f and prove a Basic Lemma. In section 3 we deal with the 
equation PA + BP = -C. In section 4 we analyze the Lyapunov equation, 
give algorithms for its solution and comment on the arithmetic complexity. 

In section 5 we deal with the equation P - A 'PA = Q and in section 6 we 
present numerical examples. 


V 
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2. The action f„_ 
— - BA 


Let A be an nxn matrix and B an mxm matrix both over the reals. Let 
R[x,y] be the ring of polynomials in two indeterminate x and y over the 
real numbers R. Let V - (<J> 2 (x), ij> 2 (y)) be the ideal in R[x,yl generated 
by <|> 2 (x) the characteristic polynomial of A, and \J> 2 (y) the characteristic 
polynomial of B. Elements of the quotient ring R[x,y]/¥ are cosets 
denoted by ¥ + a(x,y). 

Define the action f :R[x,y] x MW -► MW in the following manner: 

Bn 


f BA (h(X,y)#M) 


I h,.B^*M*A k 
jk ^ 


where h(x,y) = ][ h x k y^, is an element in R[x,y] and M an element in MW. 


jk 


i) 

ix) 


It can be shown [2] that f has the following properties: 

BA 

f BA (u,M) * M where u is a real number 
f BA {g{x ' y)+h(x,yi ' M) " f RA ( 9<x,y) »M) + f nA (h(x,y) ,M)) 


BA 


iii) f BA (g(x*y) *h(x f ;) ,M) - f^ (g(x,y) ,f BA (h(x,y) ,M) ) 


f BA (h(X,y ) ' f BA (g(X,y > ,M) ) 


iv) Let g(x,y)mod¥ denote the polynomial of minimal degree in H* + g(x,y) 
(which cam be found by first dividing g(x,y) by 4> 2 (x) # obtaining 
the remainder R x (x,y) and in turn dividing R x <x,y) by \j> 2 (y) and 
picking its remainder) . 

Then: f^tg (x # y) ,M) - f BA (g(x # y)mod'i', M) 

v) f BA <9(x ' y)# M+N> “ * f BA (g(x,y) ' N) 

for all g(x,y) in R{x,yl and M,N in MW. 


The definition of f allows for the interpretation of MW as an 

BA 
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Rtx,y3A- nodule. 

Basic Lemma . The set MN of m*n matrices with real entries is a module 
over the quotient ring R[x,y]/1'. 

Proof of Lemma: The set of m*n matrices under addition is an abelian 

group. Define multiplication (*) of cosets ¥ + h(x,y) and m*n matrices M 
in the following manner: 

(¥ + h(x,y))*M » f BA (h(x,y)mod4', M). 

The multiplication is well defined and satisfies the properties: 

1) (¥ + h(x,y))*(M + N) * (f + h(x,y))*M + (4* + h(x,y))*N 

2) (4* + h(x,y) ) *t + g(x,y) ) *M] = [O' + h(x,y))*W + g(x,u))]*M 

3) [(4* h(x,y) ) + (4 f + g(x,y))]*M « (4' + h(x,y))*M + (4' + g(x,y))*M 

4) (4' + 1)*M » M 

for all M,N in MN and all V + h(x,y),4 f + g(x,y) in R[x,y]/4' , with 4'+l 
being the multiplicative identity in Rlxjyl/f'. 

Property v) of the action guarantees 1. Property iii) ensures the 
validity of 2. Property ii) makes certain that 3 holds. Property i) 
ensures the correctness of 4. 

3. The equation PA + BP » -C 

The Basic Lemma provides the groundwork for the construction of a 
method for obtaining the solution P of the equation 

PA + BP * -C (3.1) 

whenever a unique solution does exist. 

Equation (3.1) can be written as 
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f (x+y,P) - PA + BP - -C. 

BA 

Suppose that there exists a coset 4* + ^(x^y) such that 

(Y + q^tx^y) ) • (4* + (x+y) ) ■ (V + u) (3.2) 

where f * u is a coset which contains the real number u >< 0. Let 
q^tXfyJmod'F be the polynomial of minimal degree in ¥ + q^(x,y). We 
then have i 

(4* + (x+y) ) *P - PA + BP - -C 

t (4* + q u (x,y))* (V + (x+y) ) 1*P ■ (4* + qjx.y) ) * (-C) 

(Y + u)*P » (H' + q u (x,y))*(-C) 
uP « (4* + q u (x,y) )*(-C) 

P • U <) 

The idea therefore is to ensure that for 4* + (x+y) condition (3.2) 
holds and to then construct such a polynomial q u (x,y)mod4'. 

Proposition 1 . The coset 4* + (x+y) contains the polynomial x+y. There 
exists a coset 4* + q u (x,y) for which we have 

(4* + q u (x,y))-(4' + (x+y)) « 4' + u (3.3) 

where 4* + u is a coset containing a real number u ft 0 if and only if 

^ 0 where X^, i £ i <_ n are the eigenvalues of 4> 2 (x) c det(lx-A) 

and y , 1 <_ j < m are the eigenvalues of $_(y) ■ det(Iy-B). 

3 — i 

Proof of Proposition; We prove this Proposition by first showing that 
X t + yj ^ ° for all i,j iff ^(x) ■ vp 2 (-x) and 4> 2 (x) are relatively 
prime. Assume that ^(x) and $ 2 (x) are relatively prime. Suppose then 
that there exist X i# such that + yj - 0. This means that - -y^ 
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which implies that ^(x) and $ 2 (x) have at least one root in common. 

This in turn implies that <J> 2 (x) have a non-trivial common divisor 

which is a contradiction. Assume on the other hand that X^ + p^ + 0 
for all i,j. Suppose then that there exists a k(x) of degree greater than 
or equal to one such that k(x)|if^(x) k(x)j4> 2 (x). This would imply that 

(x) and $ 2 <x) have at least one root in common which contradicts our 
assumption. 

It can be shown [2] that 
x+y|$ 2 (x)ij> 2 (y) - 0 1 (y)^ 1 (x). 

<Mx)\|> (y) - <J> (y)^> (x) 

Let p w< x »y> * (3,4) 

We now prove the Proposition. 

Assume that X^ + p^ t 0. We then have that 4> 2 (x) *nd ^(x) are 

relatively prime# which implies that there exist polynomials X^tx) , p g (x ) , 

X'(x)#M'(x) such that 
e e 

\.(x)*fi.(x) + p (x) <J>_ (x) « e 

61 e ^ ( 3.5) 

XMx)ij; 2 (x) + p^txj^tx) ■ e 

for some element e ft 0 in R. 

Let q u (x,y) - X e (x)Ug(y)P^(x,y) . 

Since 

(x+y) *q u (x,y) « X e <x)p^(y)P^(x,y) 

= X 0 (x) p^ (y) (j> 2 (x) ^ 2 (y) + eXMy)i^ 2 (y) 

+ ep^ (x) $ 2 (x) - p e (x)X^(y)<$> 2 (x)t^ 2 (y) ~ e 2 . 


2 

we must have (u = -e ) 
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> 


(f + (x4y))*(f + q u (x f y)) ■ ¥ + u. 

Assume on the other hand that there exists a coset ¥ + q^(x,y) such 
that (¥ + q u (x,y) ) * (¥ + (x+y)) ■ ¥ + u where ¥ + u contains the real 
number u f< 0. Show that X^ + f 0 for all i,j. 

We have that 

q u (x,y) • (x+y) - a(x,y)$ 2 <x) + b(x,y)ij> 2 (y) + u. (3.6) 

Suppose that there exist i ■ i' and j ■ j' such that 



Evaluating (3.6) at x ■ X^, and y ■ y^ f we have that 
0 ■ u 

which is a contradiction. This completes the proof of Proposition 1. 


As can be seen from the proof of Proposition 1 the polynomial q^fx.y) 
can be constructed and this prescribes an algorithm for the solution of 
equation (3.1) . 

Algorithm for solving the Linear matrix equation PA + BP * -C. 

Al) Obtain $ 2 (x) , ; (x) the characteristic polynomials of matrices 

A and B respectively. 

$ 2 (x)i|> ( y) - <My)iMx) 

A2) Set P. . - ~ . 

{(4 X + y 

A3) Using the Extended Euclidean algorithm or an equivalent method 
obtain the polynomials X e (x), uMx) and e. 


A4) Find q u (x,y) «X e (x)vi^(y) P^U»y> • 
A5) Form P «* f (q (x,y)mod¥, -C) . 

U DA U 

A6) Set P » ■— P , u » -e^. 


u u’ 
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4, The Lyapunov equation A'P + PA 


Suppose that B ■ A‘, C • Q ■ Q' with A a stability matrix (one 
which has the real parts of its eigenvalues in the left half complex 
plane). This is the special case of equation (3.1) known as the Lyapunov 
equation. Because of xts importance we study it separately. 

In this case where B • A' let us denote the action f^ by f^. Let 
* ($ 2 (x), $ 2 (y)) be the ideal in R[x,y) generated by $ 2 (x) an<S $ 2 W 
where as previously $ 2 (x) is the characteristic polynomial of A. We 
denote by t + q(x,y) the cosets in R[x,y}/$. We then have the following 
corollary to Proposition 1. 


Corollary 1 : The coset $ + (x+y) contains the polynomial x+y. There 

exists a coset $ + q^(x,y) for which 

($ + q^tx.y))*^ + (x+y)) ■ $ + u 

where $ + u is a coset containing the real number u ft 0, if and only if 
X^ + X^ ft 0 for 1 £ i» j £ n where X^, 1 < i < n are the eigenvalues of 
$ 2 U) . 

In this case we have 


$l(x) » * 2 <-x) 


(4.1) 


T e ( x > <^1 < x ) ♦ X e (x)$ 2 (x) ■ e 


for a real number e ft 0 


P^(x,y) 


$ 2 (x)* 2 (y) - 
x + y 


(4.2) 


q u (x,y) * T e (x)T e (y)P^(x,y). 


(4.3) 
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Algorithm for solving the Lyapunov equation A'P + PA ■ -Q. 

R^) Obtain $ 2 (x) the characteristic polynomial of A 

4> 2 (x)4> 2 <y) - <My)$ <x) 
ty Set P$(x,y) - . 

Using the Extended Euclidean algorithm or an equivalent 
method obtain T (x) and e. 

R 4 ) Form q u (x,y) - T ft (x) t e (y) P^ (x,y) . 

R & ) Find P u - f A (q u (x # y)mod4', -Q) . 

R fi ) Set P - J • P u , u - -e 2 . 

Computer Implementation 

Since we are interested in an exact computer solution we restrict 
the field of interest to that of the rational numbers F. The algorithm 
is fully implementable, using the remarkable facilities provided by the 
computer programming system MACSYMA available at M.l.T. MACSYMA is a 
large computer programming system used for performing symbolic as well 
as numerical computations. 

Three versions of the algorithm have been constructed and p- rammed 
on MACSYMA. They are the Rational algorithm, the Integer Algor .thm and 
the Modular Algorithm having names indicative of the mode in which arith- 
metic operations are carried out. 

The Rational Algorithm 


It consists of carrying out steps through R^ in rational arithmetic. 
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The Integer Algorithm 

Suppose that the matrices A and Q only contained integer entries. 
The polynomials $ 2 (x), P$(x) then have integer coefficients. Define S 
to be the n*n matrix. 

|~a^ Sq 0 0 0 0 ...o'! 



n n*! 

where <t> 2 (x) = a^x + a^x + ... + a n and a^ = 0 for k > n. Since 
$ 2 (x) is a stability polynomial det s > 0 II], and a^ 0 < i < n are 
positive integers. If we let e = 2detS the linear system 



H“1 

has an integer solution and there exists a polynomial T^fx) * d^x + 

d„x n ... + d with d. ■ (-1) c. which satisfies 
2 n l l 

T (x)<J>.(x) + X (x)<f> 0 (x) = e. 
el e 2 

This means that q u (x,y) as in (4.3) has integer coefficients. The polynomial 
q u (x,y)mod<I> also has integer coefficients which implies that 


P u * f A<Si (X ' y) '’ Q) 


is a matrix with integer entries. 
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The algorithm proceeds as follows. 


I l ) 

Find <f> 2 < x ) 

the characteristic polynomial 



<f> 2 (x) ^2 “ <*>! ^ Cy) 

V 

Set P^fxyy) 

x + y 

I 3 ) 

Find T (x) 
e 

and e. 

V 

Form q u ix,y) » ^(xJ^CyJP^fxjy) . 

V 

Find P « f 

[ A <q u <X ' y)mod *' 

V 

Set P « — • 
u 

_ 2 

V u » -e . 


The Modular Algorithm 

The integer algorithm paves the way for a modular approach to the 
. olution. 

Suppose p is a prime that does not divide e = 2detS. If A = (a^) 

and Q = (q. .) are matrices with integer entries let Q = (q. .modp) and 
lj P 13 

A = (a. .modp) be considered as matrices over Z . A left subscript p 
P 13 P 

on a polynomial b(x,y) written as b(x f y) denotes coefficient reduction 

P 

modulo p. Suppose that coefficient arithmetic is done modulo p. We then 
have 

<f>_(x) = det (lx - A) 

p 2 p 

- . . ' d V X W y) 

pV*' 7 ’ = rTJ 

T (x) <J>.(x) + X (x) <t> (x) = e 
p e pi p e p T 2 p 

p q u (*.y) = p T e (, Ve ( y ) pV x 'y ) - 


Let P = f. ( q (x,y)mod - Q) where all arithmetic is done modulo p 
p u A p*u p p 

and $ = ( d>_(x), <f>_ (y) ) in Z [x,y] . If P and u are obtained for a 

P P 2 p 2 p pup 
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sufficient number of .primes, the Chinese Remainder Theorem can be used 
to find P and u making it possible to obtain the solution P ** — • ? . 


u u 


The algorithm is as follows: 

M.) Obtain A, Q. 

l p p 

M„) Obtain 4> 0 (x) = det{Ix - A), 

z p z p 

d> (x) (f> (y) - <j> (x) 4> (y) 

m 3 > Set P V x ' y) - E — E — — 6 


M J Obtain T (x) , e . 

4 pep 

M c ) Set a (x,y) = T (x) T (y) P A (x,y). 


P vi 


p e p e 


M.) Obtain P. = f { q (x,y)mod $, - Q). 


P u 


A p^u 


M 7 ) Repeat steps 1-6 for a sufficient number of primes and 

2 

and using the Chinese Remainder Theorem find P^ and u ■ -e . 


M Q ) Set P = - • P . 

8' u u 

Since considerable coefficient growth takes place in intermediate 
computations of the Integer algorithm a lot of storage is being used up. 
In such cases it is advantageous to use the Modular Algorithm. 


Arithmetic Complexity of the Integer Algorithm 

We are concerned with the number of integer operations (addition, 
subtraction, multiplication, division) involved in running the Integer 
Algorithm when A and Q are n*n matrices, using classical operations. 

Step 1^: There are several methods for obtaining the Characteristic 

polynomial <f> 2 (x) of a stability matrix. Evaluating <J> 2 (x) at n distinct 

4 

points and then solving for the coefficients requires 0(n ) operations. 
If n is in the range n £ 20 evaluating <J> 2 (x) at x=l where <f> 2 (l) = A, 


r 


r 


I 


I 
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X = flog 10 A"|and then at x = 10^ allows one to "read off" the coefficients 

3 

of <\> 2 (x) from a large integer. This procedure requires only 0(n ) 
operations. 

2 

Step 12* This step can be done in 0(n ) operations. 

Step 13; Solving a linear set of n equations simultaneously is an 0(n 3 ) 
operation. 

Step I4: Performing the multiplication as T e (x) tt^ty) 3 requires 

0(n 3 ) operations. 

Step I5: Obtaining q u {x,y)mod$ involves two polynomial divisions and 

can be done in 0(n 3 ) operations. To form f» (q (x,y)mod<5>, -Q) we use 

A U 
4 

0(n ) operations. In the event that the matrix 0 is a product of vectors 

Q = C’C' this calculation can be done in 0(n 3 ) operations. 

2 

Step I& : It can be done in 0(n ) operations. 

4 

It can therefore be seen that the overall calculation requires 0{n ) 
operations in general and 0(n 3 ) operations in the special cases mentioned. 

storage requirements are much harder to determine since the imple- 
mentation is on a variable length word computer. 

5 . The equation P - A* PA * Q 

The Basic Lemma provides the groundwork for the construction of a method 
for obtaining the solution P of the equation 

P - A’ PA - Q (5.1) 

whenever a unique solution exists. 

Equation (5.1) can be written as 
f^d - xy, P) = P - A' PA * Q 


(5.2) 
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where as in the Lyapunov equation, we are interested in the case when 
B ■ A' and we denote f by f ft . Let $ * (4>2<x) , <p^(y)) with <J> 2 (x) * 
det(Ix -A). We then have the following Proposition. 

Proposition 2. The coset $ + (1 - xy) contains the polynomial 1-xy. 
There exists a coset $ + q u (x,y) for which 
(<!> + q u (x,y))d> + (1-xy)) = $+u 


where $+u is a coset containing the real number u^O, if and only if 
1— A^A . ft 0 for 1 < i,j < n where A i 1 i n are the eigenvalues of 


<P 2 (x) . 

If we have that 

, . , n , n-1 , 

(x) = a x + a -X + ... + a A 
T 2 n n-1 0 


then let 

a , x n, . -1. n , n-1 

0 3 (x) = x <t> 2 (x ) = a Q x + a^x + ... + a n - 

We can see that the roots of <t> 3 (x) are 1 < i < n. 

i 

Proof of Proposition 2; We will first show that 4> 2 (x) 3,1(5 <t> 3 (x) are rela- 
tively prime if and only if 1 - A. A. / 0. Assume that 4>, (x) ,<p 0 (x) are 
relatively prime. Suppose then that there exist A^,A^ such that 1 - A^A^ = 0. 
This means that A^ = which implies that (x) and 4> 2 <x) have at least 
one root in common. This implies in turn that (f> 2 (x) ,<J> 3 (x) have a non- 
trivial common divisor which is a contradiction. Assume on the other hand 
that 1 - A^Aj ^ 0 for all i,j. Suppose that there exists a k(x) of degree 
greater than or equal to one such that k(x)|<f> 2 (x) and k (x) | <J> 3 (x) . This 
implies that 4> 2 (x), <j> 3 (x) have at least one root in common which contradicts 
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our assumption. 

It can be shown that 
1 - xy|0 2 (x)4» 2 (y) - ^ 3 (x)<|» 3 (y) . 

Let 

4> 2 < x >$ 2 <y> * 4> 3 (x)$ 3 (y) 

V*' y) • 

We now prove the Proposition. 

Assume that 1 - X^X^ l* 0. We have that 4> 2 (x), 4*3 (x) 

prime which implies that there exist polynomials T e {x), 

T (x)<Mx) + X (x) 4> 0 <x) * e 
e J e 2 . 

for some element e ft 0 in R. 

Let q u (x,y) = T e (x)T e (y)P^(x,y). 

Since 

q u (x,y) (1-xy) =* T e (x) T e (y) P^ (x,y) (1-xy) 

= T (x) T (y)4> (x)4>, (y) 
e e z & 

+ eX e (y)4> 2 (y) + eX e (x)<J> 2 (x) 

- X (x)X (y) 4> 0 (x) (j). (y) - e 2 
e e z z 

2 

we must have (u = -e ) 

(<t> + (l-xy))*(4> + q u (x,y)) = $+u. 

Assume on the other hand that there exists a coset $ 

that 

(4> + q u (x,y))($ + (1-xy)) = $ + u 
where 4>+u contains the real number u 7* 0. Show that 1 - X 


are relatively 
(x) such that 


+ q u ( x , y) such 

. X . ^ 0 for all i, j . 
1 j J 


We have that 
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q^Xfy) • (1-xy) - a(x,y)<)> 2 (x) + b(x,y)<J> 2 (y) + u (5.3) 

* 

Suppose that there exists i « i' and j ■ j' i.uch that 1 - X^ f Xj, - 0. 
Evaluating (5.3) at x * X^,, y ■ X^, we have 
0 - u 

which is a contradiction. This completes the proof of Proposition 2. 

As can be seen from the proof of Proposition 2, the polynomial q u (x,y) 
can be constructed and this prescribes an algorithm for the solution of 
equation (5.1). 

Algorithm for solving the linear matrix equation P - A' PA ■ Q. 

B^) Obtain $ 2 (x) = det(Ix - A). 

<Mx)<j> (y) - <J> (x)4> (y) 

B ) Set P A <x»y) - — ~ - - . 

2 <P 1-xy 

B^) Using the Extended Euclidean Algorithm or an equivalent 

method obtain T e (x) * e. 

B 4 ) Form q u (x,y) « T e (x> T e (y)P^(x,y) . 

B 5 ) Form P u =» f ft (q u (x,y) ,Q) . 

B.) Set P = - • P . 

6 u u 

6. Numerical Examples 
We wish to compute 
00 

G ** f x' (t) *Q*x(t) dt 
J 0 

where x(t) is a solution to 

x(t) = Ax(t) x(0) * c. (*) 


The system modelled by (*) is of the form 



where the number of blocks is finite. 

Example Is The number of blocks is 5 with C-l, K«l, M-10000. Listed are 
the corresponding A matrix, the Q matrix and the solution P to the 
equation PA + A'P » Q. 

Example 2: The number of blocks is 2. Listed are the corresponding A 

matrix in parametric form, the Q matrix and the parametric solution P of 
the equation PA + A'P ■ Q. The parametric solution P is valid only for 


appropriate values 

of E, 

M, Z, 

(Z-O. 






Example 1: 
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5000 
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0000000000 
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